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Summary 

In an effort to support the design of zirconium diboride 
composite leading edges for hypersonic vehicles, a finite 
element model (FEM) of a prototype leading edge was 
created and finite element analysis (FEA) was employed 
to assess its thermal and structural response to aero- 
thermal boundary conditions. Unidirectional material 
properties for the structural components of the leading 
edge, a continuous fiber reinforced diboride composite, 
were computed with COSTAR. These properties agree 
well with those experimentally measured. To verify the 
analytical approach taken with COSMOS/M, an inde- 
pendent FEA of one of the leading edge assembly 
components was also done with COSTAR. Good agree- 
ment was obtained between the two codes. Both showed 
that a unidirectional lay-up had the best margin of safety 
for a simple loading case. Both located the maximum 
stress in the same region and ply. The magnitudes agreed 
within 4 percent. Trajectory based aerothermal heating 
was then applied to the leading edge assembly FEM 
created with COSMOS/M to determine steady state 
temperature response, displacement, stresses, and contact 
forces due to thermal expansion and thermal strains. 
Results show that the leading edge stagnation line tem- 
perature reached 4700° F. The maximum computed failure 
index for the laminated composite components peaks at 
4.2, and is located at the bolt flange in layer 2 of the side 
bracket. The temperature gradient in the tip causes a 
compressive stress of 279 ksi along its width and 
substantial tensile stresses within its depth. 

Nomenclature 

Ch heat transfer coefficient 

Cp specific heat at constant pressure 

h enthalpy, J/kg 

q heat flux, W/m^ 

q rac l radiant heat flux, W/m^ 


"Eloret Thermosciences Institute, Palo Alto, California. 
*MSNW, Inc., San Marcos, California. 


Sj *2 positive in plane shear strength in the 

material x-y plane 

S |2 negative in plane shear strength in the 

material x-y plane 

T temperature, K 

Q 

Xj compressive strength in the material 

longitudinal direction 
T 

Xj tensile strength in the material longitudinal 

direction 
C 

X 2 compressive strength in the material 

transverse direction 
T 

X 2 tensile strength in the material transverse 

direction 

G\ first principal stress 

02 second principal stress 

%\2 shear stress in x-y plane 

Subscripts 
r recovery 

s symmetric 

w wall 

0 stagnation conditions 

Introduction and Approach 

The present space shuttle design is limited in performance 
by the relatively low lift-to-drag (L/D) ratio of its airfoils. 
Its leading edge and nose cap are constrained to blunt 
geometries so that operational temperatures will not cause 
excessive material degradation or ablation in multiple 
reuse. Materials with higher temperature capability and 
greater ablation resistance will allow spacecraft designers 
to develop airfoils with sharper leading edge radii, higher 
L/D ratios, and better performance. 

Ultra-high temperature ceramic (UHTC) materials are 
being investigated for leading edge applications in 
advanced airframes at Ames Research Center (ref. 1). 
These materials have been evaluated in the Ames arc -jet 
test facilities, which can reproduce flow conditions 



representative of hypersonic flight. A prototype leading 
edge assembly is shown in figure 1. It was designed to 
demonstrate, quickly and simply, the functional perfor- 
mance of a sharp leading edge using existing UHTC 
materials. It is specifically designed for testing in an arc 
jet, although its size and configuration may be directly 
applicable to hypersonic vehicles. An exploded view of 
the assembly is shown in figure 2. The illustration shows 
two UHTC tip inserts, components 3 and 4. In the arc -jet 
test, or on a hypersonic vehicle, it is highly desirable to 
secure the tip inserts with a hot structure rather than an 
actively cooled one. A passively cooled attachment will 
prevent severe thermal gradients in the tip material, 
significantly reduce thermal protection system (TPS) 
weight and complexity, and increase its reliability and 
durability. Therefore, the four brackets, components 1, 2, 
5, and 6, are fabricated from a continuous fiber reinforced 
UHTC composite. Although still in development, this 
material is designed to be a load bearing hot structural 
material which is both thermochemically and thermo- 
mechanically compatible with the tip material. The arc -jet 
facility attachment, component 7, is a metallic structural 
attachment for the arc -jet facility sting. 

Initially, it was assumed that the top and bottom brackets, 
components 5 and 6 in figure 2, were the most vulnerable 
of all the components in the assembly since expansion 
of the tip inserts will force the brackets outward and, 
because of the bolted constraints, will create bending 
stresses at the root. Accordingly, the top bracket was 
modeled with the finite element method using COSTAR 
and COSMOS/M finite element analysis (FEA) codes. 
COSTAR is particularly well suited for the analysis of 
structures which incorporate composite materials but is 
not versatile enough to handle large multicomponent 
assemblies. Additionally, FEA using COSTAR with 
similar ceramic matrix composites (CMC) has been 
demonstrated (ref. 2). The COSTAR analysis was 
completed under contract by MSNW, Inc., San Marcos, 
California, and is reported in reference 3. The reason for 
the duplicate analysis was to lend credibility to our 
COSMOS/M analysis with composite materials, which 
had not been done previously. 

The fiber volume fraction and fiber orientations of the 
bracket material were important considerations to ensure 
that the brackets had adequate structural integrity to 
survive the severe heating environment. Several 
orientations and fiber loading were evaluated using 
COSTAR, and directional material properties were 
calculated. Iterations on these calculations were made 
until reasonable agreement was obtained with experi- 
mentally measured values. These properties were then 
used in both COSTAR and COSMOS/M for FEA. The 
results of these two independent analyses performed on 


the top brackets were compared and noted before 
proceeding with FEA of the leading edge assembly. 

For each structural analysis of layered composites 
presented in this report, stress results are presented in 
terms of Tsai-Wu failure criteria. Similar to a safety 
factor, the Tsai-Wu criteria (F.I.) are a function of 
principal stresses and ultimate strengths of the material. 
The Tsai-Wu criteria used here are defined by 
equation (1). 

F.I. = (F|CJ| + F 2 a 2 + F6t 12 ) 
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A value of 1.0 denotes the onset of failure. These criteria 
were calculated for each layer in the composite. The 
highest values tend to be in the outer two layers of the 
composite. 

Modeling and analyzing the entire leading edge assembly 
serves to pinpoint the exact areas of contact governed 
by the interaction of all components within it as it is 
thermally loaded and also identifies locations and 
magnitudes of thermally induced stresses in each of the 
components. But it is the contact loads that determine the 
structural integrity of the attachment concept. With this 
information one can assess whether the materials and 
design are compatible and the design concept is 
thermostructurally viable. This type of analytical 
approach will be instrumental in future designs of similar 
TPS. 


The assumptions regarding the boundary conditions used 
in the COSMOS/M and COSTAR analyses of the top 
bracket, region of contact and local temperature, were 
later refined in modeling the leading edge assembly. 

A trajectory based aerothermal boundary condition was 
obtained by selecting a coordinate on the tip materiafs 
aerothermal performance constraint which coincided with 
conditions obtainable in ground test facilities (arc jets) 
and with a coordinate on a projected hypersonic vehicle’s 
flight trajectory (ref. 4). From this initial condition, 
boundary conditions were calculated as a function of 
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streamline distance so that the heating distribution over 
the entire assembly could be properly treated. 

Analyses and Results 

Composite Materials and Properties 

In order to properly design the leading edge, it is 
necessary to understand the material behavior and com- 
pute accurate stresses in its components caused by the 
operating conditions. The objectives in analyzing the 
composite material behavior are to (1) develop a realistic 
micromechanical model of the composite material 
properties, (2) compute properties for composites made 
with various fiber volume fractions, (3) select an 
optimized laminate design for fabrication and testing. 

The tip inserts, components 3 and 4 in figure 2, are 
fabricated from a ZrB 2 /SiC particulate composite with 
SiC as the minor phase. Billets of this material are 
processed by unidirectional hot-pressing and the inserts 
are machined from the center of a billet. Material 
properties for ZrB 2 /SiC were obtained from the Ames 
TPSX database and were assumed to be isotropic for 
these analyses. Processing of the ZrB 2 /SiC material 
has been refined by White Materials Engineering of 
Cumberland, Rhode Island, under NASA ARC Contract 
No. NAS2- 14242 and is being marketed under the trade 
name “Zirstar.” 

The brackets, components 1,2, 5, and 6 in figure 2, are 
fabricated from a continuous fiber reinforced di boride 
particulate composite. The bracket material matrix is of 
the same composition as the tip inserts but is reinforced 
with continuous SiC fibers (Textron SCS-9a). It is 
processed in thin sheets, or lamina, which are then 
stacked in the direction of choice and consolidated by 
hot-pressing. The material has been developed by 
Advanced Ceramics Research of Tucson, Arizona, under 
NASA ARC Contract No. NAS2- 13796. It is designated 
as SCS9/RBSiC/ZrB 2 . 

Unidirectional lay-ups, 0° fiber orientation, containing 
20.7 vol.% fiber had been characterized extensively in 
both the fiber and transverse direction. The anisotropic 
properties of this material were calculated with COSTAR. 
The calculated modulus of elasticity was adjusted until 
it agreed reasonably well with that measured over a 
temperature range from room temperature to 3000°F. The 
computations resulted in effective mechanical properties 
of the matrix that were significantly less than those 
reported for the fully dense ZiTWSiC material. One 
reason for this is because the matrix is heavily micro- 
cracked due to the thermomechanical mismatch between 
the fiber and matrix. Because experimental measurements 


of the thermal properties were not available, those of 
ZrB 2 /SiC were used. Once agreement was obtained with 
the 20.7 vol.% loading, the same effective matrix was 
used in calculating properties for composites with fiber 
volume fractions of 10, 22.5, and 35 vol.%. It was found 
that the 0° modulus and tensile strength are linearly 
dependent on fiber volume fraction, so the composite 
with 35 vol.% loading exhibits the highest strength. The 
engineering properties of this material are given in 
table 1. Because the 35 vol.% loading has the greatest 
strength, it was chosen for the top bracket bending 
analysis. 

Top Bracket Bending 

A finite element model of the leading edge top bracket 
was developed using COSTAR and COSMOS/M to 
analyze the performance of 35 vol.% SCS9/RBSiC/ZrB2 
laminates with orientations of 0°, 0°/90°, and 
0°/±45°/90°. The models were clamped at the bolt 
and loaded by a tip deflection that was representative 
of the thermal growth of the Zirstar leading edge test 
specimen. All analyses were performed at a uniform 
temperature of 3000°F. These finite element models were 
verified by comparison with beam theory calculations. 

COSTAR results include the distribution of plate forces 
and moments within the bracket and ply-level stresses 
within the composite laminate. Results of the stress 
analysis showed that the bracket is loaded primarily by a 
bending moment, which generates longitudinal bending 
stresses in the bracket. The comparison of the results for 
the various orientations showed that the unidirectional 
design (i.e., 0°) exhibited the largest margin of safety of 
all the laminates that were analyzed. The results also 
showed that the performance of the 0/90 design was 
optimized by concentrating the 0° plies at the outside 
surface of the bracket and situating the 90° plies near the 
neutral axis. Based on these results, the 0° laminate was 
selected as the most attractive material design and is 
recommended for the bracket fabrication. 

The COSMOS/M finite element model was created using 
the same top bracket geometry and material properties 
as the COSTAR model (see fig. 3 and table 1). The 
COSMOS model consisted of 890 quadrilateral composite 
shell elements as did the COSTAR model. The grid 
included half of the width of the bracket with symmetry 
conditions (UY = RX = RZ = 0) applied along the 
centerline. Figure 4 shows the COSMOS/M model mesh. 
The element coordinate system definition used for stress 
computation is also described in figure 4. The element 
coordinate system was chosen such that it matches the 
coordinate system of the COSTAR model and the 
coordinate system defined in the geometry description 


3 



(fig. 3). The element coordinate system can be visualized 
as a weave of fibers such that the longitudinal fibers 
(90° orientation or x direction) always follow the curved 
profile along the length of the bracket. Nodes at the bolt 
location were constrained in all three translational degrees 
of freedom (UX = UY = UZ = 0). 

The results for the 0/90 configuration for the outer two 
plies are shown in the form of Tsai-Wu failure criteria in 
figures 5-8. Both codes give similar results for stress 
magnitude and distribution. The highest stress in both 
models occurs in ply 2 at the bolt interface (see figs. 5 
and 8) where the failure index reaches 2. 15 for COSTAR 
and 2.23 for COSMOS. These results suggest that failure 
will occur in the bolt region of the top bracket when the 
Zirstar tip expands and bends the bracket. It should be 
noted that in both models the maximum stresses are 
unrealistically high due to the constraints imposed to 
simulate bolting the bracket (constrained in all degrees of 
freedom). For comparison, the stress for F.I. maximum is 
plotted in figure 9 for each ply through the thickness. 

Both codes produce nearly identical results from the 
outermost ply to the neutral axis for plies 1 and 2. Both 
codes also produced erratic results after the neutral axis 
for the 90° plies. The reason for this behavior is not 
known at this time. Broader contact forces and friction 
were not considered in this study but would represent a 
more realistic treatment of the bolt interface for fastening 
the structure. A detailed, more complex analysis of the 
bolt region deserves further study. 

To refine and validate some of the assumptions used in 
the COSMOS/M and COSTAR top bracket models, a 
bi-planar symmetric 3-D thermostructural model of the 
entire test fixture with trajectory based aerothermal 
heating was created with COSMOS/M. Benefits of 
modeling the entire assembly include a more refined 
temperature distribution that takes into account realistic 
thermal loading, edge effects due to reradiation, and a 
better understanding of the component contact locations 
and contact forces. Although it was found from the top 
bracket bending analysis that concentrating the 0° plies at 
the outside surface of the bracket and situating the 90° 
plies near the neutral axis is preferable for reducing stress 
in the top bracket, a 0°/90° symmetrical lay-up was used 
to formulate the full assembly analysis. This is because 
the model being tested in the arc -jet facility was fabri- 
cated with a 0°/90° symmetrical lay-up before the 
analysis was complete. A major reason for performing a 
3-D assembly analysis was to predict the survivability of 
the leading edge model. However, future prototypes (and 
analyses) will utilize a 0°/0° lay-up for the outer plies. 


Thermal and Structural Analysis of the Leading 
Edge Assembly 

Before the structural response of the leading edge 
assembly can be determined, an accurate steady state 
temperature distribution must be established for the 
condition of interest. To accomplish this, a finite element 
model (FEM) of the leading edge assembly was created, 
and heat flux was applied to the exposed surfaces corre- 
sponding to a point on an aerothermal performance 
constraint. Heat flux was applied to the model by pre- 
scribing a convection coefficient to each element along 
the streamline from the stagnation line to the sting mount. 

Thermal and structural FEM- To create the FEM used 
for both the thermal and structural analyses, symmetrical 
representations of the top and side brackets and tip inserts 
were meshed with appropriate elements. Geometrical 
and load symmetry permitted the use of symmetrical 
constraints about the xy and zx planes to make efficient 
use of the element density. The mesh assembly is shown 
in figure 10. A total of 5000 elements were used. The 
tip insert was modeled with eight node isotropic solid 
elements. An increased element density was used in the 
tip region where a high thermal gradient is expected. 

The top and side brackets were modeled with eight node 
orthotropic SHELL4L elements. The SHELL4L elements 
were assigned both thermal and structural orthotropic 
material properties with temperature dependence. A 
20 layer [0/90] s composite architecture was defined in 
the SHELL4L real constant set. Material properties were 
computed for this architecture in the same fashion as that 
done for the top bracket analysis of this report. Gap 
elements were used between components to simulate a 
finite gap on the order of 0.001 inch. As the components 
expand when heated, and the gap distance becomes zero 
within a specified tolerance, the components will be 
capable of reacting together structurally. Unfortunately, 
gap elements available in COSMOS/M do not have the 
capability to conduct heat across the gap even if the gap 
distance becomes zero. In reality, heat conduction will 
occur when the gap distance approaches zero. To work 
around this problem, very soft (E x = IE-5 psi) TRUSS3D 
elements with large cross-sectional areas (Acs = 1 in.) and 
a high thermal conductivity (k x = 1E5 Btu/h-ft-°F), were 
positioned coincident with the gap elements to simulate 
heat conduction between them. The use of the soft 
conductor bars imposes a constant contact condition 
without altering structural response. 

Boundary conditions- The aerothermal heating 
conditions were obtained by selecting a coordinate on the 
tip material’s aerothermal performance constraint which 
coincides with conditions obtainable in ground test 
facilities (arc jets) and with a coordinate on a projected 
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vehicle’s flight trajectory (ref. 4). Aerothermal perfor- 
mance constraints are computed by taking into account 
all the various modes of heat dissipation available to a 
material including reradiation, convective blocking, 
catalytic effects, and multidimensional conduction 
(ref. 5). The computation is done assuming steady state 
heat transfer and results in a material based performance 
map in altitude-velocity space. 

In order to accurately calculate steady state aerothermal 
performance constraints for UHTC sharp leading edges, 
it is necessary to analyze the thermal conduction of heat 
away from the stagnation region. Rather than developing 
custom thermal analysis software, it is advantageous 
to use commercially available FEA software tools 
(COSMOS/M) which are capable of thermal/structural 
analysis. Because these commercial FEA tools have been 
developed for application to a wide variety of generic heat 
transfer problems, there are standard procedures for 
constructing the model geometry, specifying the element 
type, and assigning material properties. However, because 
of the unique characteristics of heat transfer to a nose tip 
(3-D) or wing (2-D) leading edge in hypervelocity flight, 
there are no standard procedures for assigning the thermal 
boundary conditions to the FEA model. It is necessary to 
develop a custom procedure for assigning the appropriate 
aerothermodynamic heating boundary condition. 

Aerothermodynamic heating of a leading edge in 
hypervelocity flight has a maximum where the fluid is 
stagnated, then decreases rapidly as the fluid accelerates 
downstream. The dependence of this heating distribution 
on freestream velocity, geometry, fluid properties, and 
wall temperature was first analyzed by Lee (ref. 6) for 
simple 3-D and 2-D leading edges. More generalized 
expressions which are useful for examining complex 
geometries were developed by Marvin (ref. 7). Both of 
these approaches provide estimates of normalized heat 
flux as a function of surface distance from the stagnation 
point. For a given velocity, geometry, and fluid this 
distribution simplifies to 

q w = Ch(h r - h w ) 

4w,0 ^h,o(^r “ ^w)q 

Using h = CpT 

q w _ ChCp(T r — T w ) 
q W ,0 Ch,oCp,o(Tr - T w )q 

If Cp and T r are constant, then aerothermodynamic 
heating of a leading edge is given by 

q w _ l^r - T w ) 

qw,0 Ch,o(Tf - T w )Q 


where and T w are functions of surface distance from 
the stagnation point. Typically, the standard procedure in 
thermal FEA uses up to four types of boundary conditions 
(BC) for specifying the thermal load on the model 

1 . Temperature (T w ) 

2. Convection (Ch, T r ) 

3. Heat Flux (q w ) 

4. Radiation (q ra d) 

With these four alternatives, there are two approaches to 
specifying aerothermodynamic heating for thermal FEA. 
In one approach, using BC-2, the heating can be 
expressed from equation (4) by specifying and T w 
as a function of distance along with the constants q w> o 
and T r . In the other approach, using BC-3, the heating 
is expressed from equation (4) by specifying q w as a 
function of distance along with the constant q Wt o- In 
addition to either approach, BC-2 or BC-3, reradiation 
must be specified by assigning the appropriate surface 
emissivity using BC-4. 

For the analysis discussed in this study the aero- 
thermodynamic heating q w /q Wi o is specified using the 
aerothermal performance constraint code PERFORM 
(ref. 5). Because of the nonlinear coupling between the 
convective heat transfer from the fluid and the surface 
temperature of the material due to the reradiation, several 
iterations are required between PERFORM and 
COSMOS/M. Convergence on surface temperature 
usually occurs in less than four iterations when using the 
BC-2 approach. Convergence is much worse when using 
the BC-3 approach, and does not occur under some 
conditions. 

The applicable aerothermal performance curve is shown 
in figure 1 1(a) along with the position of the selected 
coordinate, an altitude of 197 kft and a velocity of 
29 kft/s. From this coordinate, heat transfer coefficients 
were computed as a function of streamline distance as 
shown in figures 1 1(b) and 1 1(c). These “hot wall” 
coefficients were input into COSMOS/M as the boundary 
conditions representing aerodynamic heating over the 
entire assembly. 

Other necessary thermal boundary conditions include 
those at the back wall or bolt interface. A constant 
temperature of 300°F was selected here for several 
reasons. The heat transfer coefficients computed for this 
case are nearly the same for an adiabatic back wall case; 
in both flight and arc -jet testing the assembly is attached 
to a structure which acts as a large heat sink, and 300°F 
is a nominal use temperature for many aerospace struc- 
tures (i.e., carbon/epoxy composites or metallic alloys). 

An emittance of 0.7 rather than 0.6 was used in the 
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COSMOS/M analysis for all components because this is 
the materials database value for the unoxidized surface 
at the temperature of interest. An emittance of 0.6 is 
characteristic of the oxidized surface. Using 0.7 results 
in a lowering of the stagnation line temperature by 
approximately 300°F from its design limit. The ambient 
temperature for the COSMOS/M analysis was assumed 
to be 500°F. 

The most critical structural boundary condition is the 
treatment of the assembly connections. The side and top 
brackets were considered to be securely bolted to the 
sting by constraining nodes associated with the fastener 
boundaries. It is believed to be a conservative approach 
since some slipping will occur at the bearing surface 
between the bolt head and the sting. The temperature 
profile generated from the thermal analysis and associated 
thermal strains is the sole source of the loads applied to 
the components of the assembly. The thermal strains are 
induced by temperature differential and mechanical 
constraints that restrict thermal growth. 

A script utilizing the COSMOS/M command language 
was developed to map discretized points from a curve fit 
of s/r versus heat transfer coefficient data onto the leading 
edge assembly. The script can be applied to any FEA 
model generated with COSMOS/M by substituting the 
geometry profile equations, reassigning the element 
step size and identifying the new stagnation location. 
(Documentation is embedded in the script which can be 
found in the Appendix.) 

Thermal-Figure 12 shows the steady state temperature 
distribution at the prescribed aerothermal condition 
looking from the top-side view. The peak temperature is 
at the stagnation line of the tip inserts and has a value of 
4670°F. The temperature contour lines drifting forward at 
the model edge are an effect of reradiation conditions 
imposed at the outward facing surfaces of the top bracket 
and conduction in the top and bottom brackets, which are 
twice as thick as the side brackets. The other side ol the 
bracket is not exposed to air, but is in thermal contact 
with the adjacent leading edge segment. 

Structural- Thermal strains develop in components of 
the assembly from thermal gradients and from mechani- 
cally restricting thermal growth. For the top and side 
brackets, mechanically constraining the thermal growth is 
the primary contribution to the stresses developed. The 
bulk of the stress in the tip insert is primarily induced 
from severe thermal gradients. A major assumption in the 
preliminary top bracket analysis was that the top bracket 
is loaded solely by the thermal growth of the tip insert 
in the z direction. Modeling all components shows that 
the tip insert imparts very little load to the top bracket. 
Inspection of the displacement plot in figure 13 


(a bottom-side view) reveals that the area of contact is at 
the rear of the assembly between the side and top/bottom 
brackets. The gap elements in this region are closed and 
generate a total compressive force of 383 lbf. This force is 
listed with other contact forces in table 2. Loading in this 
region creates a prying effect that causes the tip of the top 
bracket to separate slightly from the tip insert. The gap 
created from this separation is predicted to be on the order 
of 5 mils. 

The load built up in the contact region reflects the rigid 
constraints imposed to simulate the bolted connection. As 
a result, stresses in the bolt bearing area are extreme. A 
Tsai-Wu failure index of 4.1 was computed at the bolt 
region (identified by label L4) in layer 2 of the side 
bracket as shown in figure 14. If the bolt constraints were 
realistic, failure would occur in both layers at the bolt 
interfaces. 

Table 3 summarizes the critical stress areas in the loaded 
bracket (L1-L6). The Tsai-Wu failure indices for ply 1 
are shown in figure 14. Critical stress areas for ply 1 
include the top bracket’s bolted region (L3) and the large 
bend radius of the top bracket (LI). The stress computed 
at the targe bend radius approaches the material’s failure 
limit and is not considered to be conservative. 

Due to severe thermal gradients in the tip insert, 
considerable tensile and compressive stresses develop. 
Figures 15 and 16 show a stress contour plot for the 
z direction for layers 1 and 2, respectively. Because the 
tip inserts are considered to be made of an isotropic 
material the results do not vary by layer. Layer results 
apply only to the top and side brackets. The greatest 
stresses developed in the leading edge tip insert are 
compressive with a maximum magnitude of 297 ksi. The 
location of this stress is shown by the designation L5 in 
figure 17. The center area of the tip inserts is in tension 
with a maximum magnitude of 45 ksi. The strength of the 
zirconium diboride composite is 50 and 500 ksi in flexure 
and compression, respectively, so the tensile stresses in 
the tip insert, as shown here, could lead to failure. 

Discussion 

The static solution indicates that when the model is 
thermally loaded the contact forces developed between 
the side and top brackets, as a result of the thermal growth 
of the top brackets, cause the side brackets to separate 
from the tip insert. Upon separation, conduction can no 
longer occur. Therefore, the use of truss elements, while 
necessary to provide a conduction path across the gap 
elements, cannot fully simulate the real situation. The 
inability to capture conduction termination upon separa- 
tion is a limitation of the FEA code used. If the degree of 
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thermal conductivity could be defined for each truss 
element as a function of relative nodal positions, and the 
problem were solved using an incremental solution 
algorithm that solves the thermal and structural responses 
alternately, the effect could be better simulated. 

Conclusions and Recommendations 

Failure is likely to occur in the top and side brackets if the 
leading edge assembly is rigidly constrained as modeled. 
The FEA results for the top and side brackets are some- 
what conservative since some slipping will occur between 
the bolt and the sting flange. However, slipping is also 
undesirable since the displaced brackets will remain 
displaced upon cooldown. 

Stress in the top bracket’s large bend radius (LI) 
approaches the failure limit. It is strongly recommended 
to modify the sting attachment scheme by designing in 
more compliance to alleviate stress formed in the LI 
region. One design possibility is to restrain the top 
brackets in the z direction with a high temperature spring 
attached to the right and left top brackets instead of 
rigidly bolting them to the sting. When the side brackets 
begins to expand, the top brackets will slide outward in 
the slot formed by the two halves of the side brackets and 
will be tensioned such that, upon cooling, the top brackets 
will be forced back into their original position. 

Even though the maximum tip temperature is below the 
single use limit, failure may occur in the tip insert 
because large tensile stresses in the z direction approach 
the failure stress of the material. These stresses may be 
mitigated by shortening the dovetail of the tip insert. This 
result maybe conservative as well, however, because the 
insulating effect of the oxide layer is not taken into 
account. 
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Appendix 

C*File: s2mbeta.txt 10/95 

C*This script maps heat teansfer coefficients (as a function of S/R) 

C*to finite elements of an extruded leading edge model created with 
C*COSMOS/M. The heat transfer coefficient/S/R curves repesent 
C*conditions for vehicle re-entry at an altitude of 60 km and were 
C*determined by iteratively solving a 2d profile Heat transfer model 
C*in conjuction with ACE code. (ref. Paul Koldijiez) 

C* 

C*Tom Kowalski, Eloret 

C*Thermal protection Materials and Systems Branch (STM) 

C*NASA Ames Research Center. 

C*9/ 14/95 
C* 

C*NOTE: 0,0,0 IS AT THE RIGHT MOST RADIAL CENTER 
C* 

C*Function definitions for hcoef and s/r 
C*PARASSIGN, Al, REAL, -9.091 1925 
C*PARASSIGN, Bl, REAL, -1.004798 
C *P ARASSIGN, C 1 , REAL, 9.0405 1 2 
C*PARASSIGN, Dl, REAL, .37909244 
C*P ARASSIGN, El, REAL, -3.613106 
C*P ARASSIGN, A2, REAL, 6.4622922E-5 
C*PARASSIGN, B2, REAL, -3.56788E-5 
C*P ARASSIGN, C2, REAL, -2.3938027E-5 
C*P ARASSIGN, D2, REAL, 8.1591626E-6 
C*PARASSIGN, E2, REAL, 5.61 804 IE-6 
C*P ARASSIGN, F2, REAL, 6.7637685E-7 
C* 

FUNCDEF, HTC02(X1), REAL, EXP((-9.091 1925& 

+9.0405 1 2*X 1 -3 .6 1 3 1 06*X 1 A 2)/( 1 - 1.004798 *X 1 +& 

.37909244*X 1 A 2)> 

FUNCDEF, HTC270(X2), REAL, (6.4622922E-5)& 
-((3.56788E-5)*LOG(X2))-((2.3938027E-5)/LOG(X2))+& 

((8. 1 59 1 626E-6)*(LOG(X2)) A 2)+& 

((5.618041E-6)/(LOG(X2)) A 2)-((6.7637685E-7)*(LOG(X2)) A 3) 

C* 

PARASSIGN, LENGTH, REAL, 6. 80 
PARASSIGN, ALPHA, REAL, (5*PI/180) 

PARASSIGN, RAD, REAL,. 1 1 
PARASSIGN, MINESZ, REAL, .0145 

PARASSIGN, STRLINE, REAL, (LENGTH-RAD+RAD*SIN(ALPHA))/COS( ALPHA) 
PARASSIGN, STRARC, REAL, RAD*((PI/2)-ALPHA) 

PARASSIGN, STRTOT, REAL, (STRLINE+STRARC) 

PARASSIGN, TSTRSTPS, REAL, CEIL(STRTOT/MINESZ) 

PARASSIGN, STRIN, REAL, 1 
PARASSIGN, TZSTPS 1 , REAL, 4 
PARASSIGN, ZSTPSZ1, REAL, .2175 
PARASSIGN, ZSTART1, REAL, 0 
PARASSIGN, TZSTPS2, REAL, 6 
PARASSIGN, ZSTPSZ2, REAL, .33416 
PARASSIGN, ZSTART2, REAL, .87 
PARASSIGN, ZLOC, REAL, ZSTPSZ1 
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C*Face numbering issue 

C* Assign location along streamiline at the end of the first 
C*set of brick elements 
PARASSIGN, ebrksetl, REAL, .19053 
PARASSIGN, sbrkset2, REAL, .22563 
INITSEL, EL, 1, 0; 

INITSEL, ND, 1, 0; 

ACTSET, CS, 19; 

C* 

C*#LOOP LABEL 10 TZSTPS 1 
C*#LOOP LABEL20 TSTRSTPS 
#LOOP LABEL20 10 

PARASSIGN, STRLOC1, REAL,STRIN*MINESZ 
PARASSIGN, SR, REAL, STRLOC 1/RAD 
#IF (SR>0) && (SR<=2) 

PARASSIGN, HTC, REAL, HTC02(SR) 

PARASSIGN, TAMB, REAL, 47510 
#ELSE 

PARASSIGN, HTC, REAL, HTC270(SR) 

PARASSIGN, TAMB, REAL, 47400 
#ENDIF 

#IF (STRLOC 1>0) && (STRLOC 1<=STRARC) 

PARASSIGN, XLOC, REAL, RAD*COS(SR) 

PARASSIGN, YLOC, REAL, RAD*SIN(SR) 

#ELSEIF (STRLOC 1 >STRARC) && (STRLOC 1 <=.22 1 76) 

PARASSIGN, XLOC, REAL, STRLOC 1 *COS(ALPHA)-RAD*SIN( ALPHA) 
PARASSIGN, YLOC, REAL, RAD*COS(ALPHA)+STRLOCl*SIN( ALPHA) 

#ELSE 

PARASSIGN, XLOC, REAL, STRLOC l*COS(ALPHA)-RAD*SIN( ALPHA) 
PARASSIGN, YLOC, REAL, .O40+(RAD*COS(ALPHA)+STRLOCI*SIN(ALPHA)) 
#ENDIF 

PARASSIGN, LOCNOD, INT, NEARND(XLOCIYLOCIZLOC) 

INITSEL, EL; 

SELREF,EL,ND,LOCNOD,LOCNOD, 1 ; 

#IF (STRLOC 1 >0) && (STRLOC 1 <=ebrkset 1 ) 

C* FIRST SET OF BRICK ELEMENTS ON CURVED SECTION, FACE 4 C*EXPOSED 

CEL, 1,HTC, TAMB,2,ELMAX, 1 ,0 

#ELSEIF (STRLOC 1 >ebrkset 1 ) && (STRLOC I <=.22 1 76) 

C* SECOND SET OF BRICK ELEMENTS, FACE 3 EXPOSED 

CEL, 1,HTC,TAMB,3,ELMAX,1,0 

#ELSE 

C*SHELL ELEMENTS, FACE 5 EXPOSED 
CEL, 1, HTC, TAMB, 5, ELM AX, 1,0 
#ENDIF 
INITSEL, EL; 

PARASSIGN, STRIN, INT, STRIN+l 
#LABEL LABEL20 

C*PARASSIGN, REAL, ZLOC, ZLOC+ZSTPSZ1 
C*#LABEL LABEL 10 



Table 1(a). Computed mechanical properties for 35 vol.% SCS-9a-[RBSiC/ZrB2] unidirectional laminate 


Temp 

E 0 

E90 

v 0-90 

Go-90 

tens 

a 0 

comp 

CT 0 

tens 

ct 90 

comp 

ct 90 

X 

(°F) 

(msi) 

(msi) 


(msi) 

(ksi) 

(ksi) 

(ksi) 

(ksi) 

(ksi) 

75 

20.74 

12.24 

0.140 

5.58 

77.4 

154.7 

4.9 

4.3 

2.8 

1000 

20.42 

13.63 

0.136 

6.19 

97.0 

193.9 

6.7 

4.3 

3.9 

2000 

19.29 

12.80 

0.136 

5.81 

75.6 

151.2 

8.2 

4.3 

4.7 

3000 

18.91 

17.01 

0.136 

7.57 

54.3 

108.5 

8.6 

4.3 

5.0 


Table 1(b). Thermal properties for 35 vol.% SCS-9a-[RBSiC/ZrB2] 
unidirectional laminate 


Temp 

(°F) 

ao 

(10- 6 /°F) 

ago 

(10- 6 /°F) 

K x 

(Btu/in.-s-F x 10^) 

Ky 

(Btu/in.-s-F x 10^) 

75 

3.29 

3.71 

3.65 

2.70 

1000 

3.34 

3.71 

1.60 

1.30 

2000 

3.34 

3.71 

0.94 

0.63 

3000 

3.54 

3.73 

0.74 

0.39 


Table 2. Contact forces between top and side brackets 


Element 

F x 

(Lbf) 

Fy 

(Lbf) 

F z 

(Lbf) 

Gap dist. 
(in) 

Normal 

force 

(Lbf) 

1093 

0.00 

0.00 

3.81 

0.00 

3.81 

1128 

0.00 

-1.29 

0.00 

0.00 

1.29 

1137 

0.00 

-0.85 

0.00 

0.00 

0.85 

1242 

0.00 

0.00 

148.30 

0.00 

148.30 

1243 

0.00 

0.00 

230.60 

0.00 

230.60 

Total 

0.00 

-2.14 

382.70 
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Table 3. Locations of concern identified in thermostructural analysis of UHTC protoype leading edge assembly 
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Tip insert Tip center at peak Low . Compressive stress Temperature difference between tip and None. Safety factor: 

temperature L5 250 ksi, half of dovetail. Hot material expands, but growth 

compressive strength is restricted by neighboring cooler material 

which places hot material in compression. 









Preliminary Top Bracket Analysis 
Loaded by Thermal Growth of Tip Insert 


.X 


4 . 


SSSSSSSfSo 

LUUJUJUJLLJLUUJlijUJ 

ooooooooo 

89S8S88?! 

n a ^ o to rt i- oj 



5 5 


X 

a> 

c LU UU 

— o o 

a> g 

co 


CD O 



16 


Figure 5. COSTAR solution , Tsai-Wu failure index for ply 1 of top bracket (0/90). 
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Figure 6. COSTAR solution, Tsai-Wu failure index for ply 2 of top bracket (0/90). 



Preliminary Top Bracket Analysis 
Loaded by Thermal Growth of Tip Insert 
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Figure 7. COSMOS solution, Tsai-Wu failure index for ply 1 of top bracket (0/90). 






Preliminary Top Bracket Analysis 
Loaded by Thermal Growth of Tip Insert 
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Figure 9. Through thickness failure index comparison for 0/90 lay-up. 




SHELL4L, 20 layer orthotropic composite 
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Figure 10. Element plot of finite element model for leading edge assembly . 





Altitude: 60 km 
Back Wall Temp: 300 °F 
T_ - 47,400 °F 


Figure 1 1(b). Convective heat transfer coefficient vs. S/r at target point indicated in figure 1 1(a) (range: 0 < S/r < 2) 
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Altitude: 60 km 
Back Wall Temp: 300 °F 
T a mb : 47,400 °F I 


Figure 1 1(c). Convective heat transfer coefficient vs. S/r at target point indicated in figure 1 1(a) (range: 2 < S/r < 70) 





UHTC Leading Edge Assembly 
Subjected to Trajectory Based Aerothermal Heating 



Figure 12. Top-side view of steady state temperature distribution plot. 






UHTC Leading Edge Assembly 
Subjected to Trajectory Based Aerothermal Heating 



Figure 13. Underside view of displacement plot. 





UHTC Leading Edge Assembly 
Subjected to Trajectory Based Aerothermal Heating 



26 


Figure 14. Tsai-Wu failure index plot — ply 2. 


UHTC Leading Edge Assembly 

Subjected to T rajectory Based Aerothermal Heating 
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Figure 15. Stress plot in the z-direction—ply 1. 



UHTC Leading Edge Assembly 
Subjected to Trajectory Based Aerothermal Heating 
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Figure 16. Stress plot in the z-direction—ply 2. 
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